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We show that nuclear pairing Hamiltonian exhibits supersymmetry in the strong-couphng hmit. 
The underlying supersymmetric quantum mechanical structure explains the degeneracies between 
the energies of the A*' and A'^max ^ + 1 pair eigenstates. The supersymmetry transformations 
connecting these states are given. 
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' Supersymmetry concepts were utilized in several nuclear physics applications. Dynamical supersymmetries relate 
, the spectra of even-even nuclei, considered as states of a system of correlated fermion pairs approximated as bosons, 
and odd-even nuclei, considered as states of a system of such bosons plus unpaired fermions [l|, Q ■ Supersymmetric 
' considerations are also used in various applications of random matrix theories For example, invariance groups of 
, the second moments arising in the study of compound nucleus scattering are orthosymplectic supergroups [4] . Similar 
' considerations arise in random matrix models mimicking QCD phase transitions [5|. Finally, it was shown that the 
0^ , spherical Nilsson Hamiltonian of the nuclear shell model has a dynamical Osp (1/2) supersymmetry connecting the 

■ SU(3) symmetry and pseudo SU(3) symmetry 

In this article we would like to discuss a fourth application of supersymmetry techniques to nuclear physics, namely 

■ to the spectra of nuclear pairing Hamiltonian: 



^ = Y1 ™ - IG*! H ^n' ^j'Sj, , (1) 

^ [ jm 33' 

where the pairing interaction between time-reversed states is written in terms of the quasi-spin operators 

m>0 m>0 

If the pairing interaction is separable, i.e. Cjf = c*Cji, then in the strong coupling limit (|G| 3> £_,) the Hamiltonian 
of Eq. ID) can be written as 

Hsc---\G\S+{0)S-{0), (3) 

where we introduced the notation 

5-+ (0) = ^ c* S'+ and S'" (0) ^ Cj S*" . (4) 

j 3 

Supersymmetric quantum mechanics has been extensively studied 0, H, 0] ■ Supersymmetric quantum mechanics 
relates the spectra of the Hamiltonians of the form A and AA\ where A is an arbitrary operator. Clearly if one 
sets A = S^{Q), then the separable Hamiltonian in Eq. ([3]) and its supersymmetric partner have related spectra. 
However, in this note we would like to highlight a more subtle supersymmetry. To this end we first introduce the 
operator 

f = exp [-iirS^^^^ (5) 

where S^'^ = ^f''^- This operator exchanges 5'+(0) and 5*^(0) 

f^S^{0)f ^ S^{Q). (6) 
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Writing S'(i) = (1/2) (5*+ + 5-), it is easy to see that the operator T transforms the empty shell (i.e particle vacuum), 
denoted by |0), to the fully occupied shell, denoted by |0): 

m = io). (7) 

It is also straightforward to show that T converts states with N particle pairs into states with N hole pairs. 
To exhibit the supersymmetry of the pairing Hamiltonian we next define the operators 

B-=ft5'-(0), B+ = S'+(0)f. (8) 

Supersymmetric quantum mechanics tells us that the partner Hamiltonians Hi = B^B^ and H2 = B^ B^ have 
identical spectra except for the ground state of Hi. This ground state has zero energy (i.e., it is annihilated by B~ , 
or alternately by 5'~(0)). The states of the Hi = B+B~, l^i), and H2 = B~ B+ , {^2), are related 

1*2) - B-\^i) = f tS-(O)l^'i) = S'+(0)f t|*i). (9) 

Using Eq. ([S]), one can easily show that the two Hamiltonians Hi and H2 are actually identical and equal to the 
pairing Hamiltonian in the strong-coupling limit: 

13+13- =13- B+ = S+{Q)S-{Q), (10) 

hence the role of the supersymmetry is to connect the states |5'2) and Clearly even though these two Hamiltonians 

are identical, the corresponding eigenstates are not. For example the particle vacuum, |0), is annihilated by B~ , 
but not by B+ . Below we show that \'^i) represent the particle states and [^2) represent the hole states. Hence 
supersymmetry connects particle and hole states. 

The energy eigenvalues and eigenstates of the Hamiltonian in Eq. ^ are worked out in detail in Refs. [l3| and 
[m and sketched in the Appendix. 

Let us first consider one-pair states. The state in Eq. p6p of the Appendix is annihilated by B~ ^ hence it has no 
supersymmetric partner. Using Eq. ([9]), we can find the supersymmetric partner of the non-zero energy state given 
in Eq. We obtain 

B-5+(0)|0) = f t[5-(0), ^+(0)]|0). (11) 

Since the action of the commutator on the particle vacuum gives a number we find 

B-5+(0)|0) -f^|0) = |0), (12) 

i.e., the supersymmetric partner of the one-pair state is the completely filled state (state with A^^ax number of pairs). 
In general the supersymmetric partners of N-pair states are states with A^max — A^ + 1 pairs. To see this consider the 
state given in Eq. with (A^ — 1) hole pairs (or equivalently with (A^max — A^+1) particle pairs). Its supersymmetric 
partner should be given by 

B+ S- {z[^^)S- (z^2^'^) . . . 5~(z^„\)|0) = S+{Q)fS-{z[^^)f^fS-{zf^)f^ . . .fS- {z^j^\)f^f\0), (13) 
which, using Eqs. ^ and (O, reduces to 

^+(0)^+(zW)...5+(zjf_\)|0). (14) 

This is nothing but the state with N-pairs given in Eq. ((22)) of the Appendix. Hence, for A^ < A^niax/2, the states with 
N pairs and (A^max — A'^+ 1) pairs have the same energy as was shown in Ref. ^l7| using an entirely different approach. 
The resulting supersymmetric spectra of the nuclear pairing Hamiltonian are illustrated in Fig. 1. Clearly this 
supersymmetry is broken by the single-particle energies (i.e. the Hamiltonian in Eq. ([3]) exhibits this supersymmetry 
while the Hamiltonian in Eq. ([1]) does not). 

We showed that nuclear pairing Hamiltonian exhibits supersymmetry in the strong-coupling limit. Provided that 
A^ < A^,„ax/2, the states with A^ pairs and (A^^ax — A^ + 1) pairs are supersymmetric partners of each other. The 
underlying supersymmetric quantum mechanical structure explains i) Existence of the zero-energy states when the 
number of pairs are less than or equal to the half of the maximum allowed value, ii) Degeneracies between the energies 
of the A^ and A^max — A^ + 1 pair eigenstates. 
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FIG. 1: Spectra of nuclear pairing Hamiltonian exhibiting supersymmetry. States with A'' pairs and with (A'^max — N +1) pairs 
are supersymmetric partners of each other (A'^ < A'max/S). 
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Appendix 



Here we summarize the information about the eigenvalues and eigenstates of the Hamiltonian of Eq. Full 
details can be found in Rcfs. and One first introduces the operators 



(x) ^ y ^ S+ and S-{x)^y- f^S", 



(15) 



where x is, in general, a complex parameter. Clearly the operators given in Eq.(j4]) are these operators with a; = 0. 
For a single pair there are two kinds of states. The state 

^+(:e(i))|0) (16) 

is an eigenstate of the Hamiltonian of Eq. ([3]) with zero energy, if x^-^^^ satisfies the Bethe ansatz equation 



(17) 



There are as many states with zero energy as different solutions of Eq. p?)l . In addition the state 

^+(0)10) (18) 



is an eigenstate with energy 



(19) 



If the available orbitals are less than half full (or at most half full) these results generalize in the following manner: 
The state 



S+ix[''^)S+{x^f>)...S+ix^^'>m 



(20) 
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is an eigenstate with zero energy if the Bethe ansatz equations 

N 



V ^ V i (21) 

^i/i„i2 AN) {N) (AT) ^ ' 



are satisfied. In addition the state 



is an eigenstate with energy 



^+(0)^+(zr)---^+(zrJi)|0) (22) 



i? = -|G||^f7,|c,p-^'4) I' (23) 

fc=i ^fe 



if the following Bethe ansatz equations are satisfied: 

51 , (W) = ~JN) + XI jlV) JAM- (24) 
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When the orbitals are more than half full there are no zero energy states. The completely filled state |0) has an energy 
given by Eq. (fT9|) . In addition the state 

^-(zr)5-(zf^)...^-(zr_\)|0) (25) 

has the same energy as in Eq. with the parameters zf^ satisfying the Bethe ansatz equations of Eq. It 

should be noted that, although the pairing problem outlined above is exactly solvable, it is not easy to find solutions 
of the resulting non-linear Bethe ansatz equations. Several strategies to solve these equations are discussed in Refs. 
[H and [H. 
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